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I. INTRODUCTION 



In the statistical mechanics of fluids^ the hquid stated is a particularly fascinating one. A 
liquid is neither a gas nor a solid, but the state where correlations really play an important 
role. The pioneering work of Berni J. Alder- showed that, because of the absence of attractive 
forces, the hard-sphere fluid admits only a single fluid phase. In order to find the liquid phase 
it is sufficient to add an attractive square-well to the pair-potential of the hard-spheres. The 
resulting hard-sphere square-well fluid admits a bell-shaped gas-liquid coexistence curved 
with the critical point moving at low temperatures and high densities as the attractive well 
width diminishes. Recently N. Kern and D. Frenkel^ studied, through computer experiments, 
a new fluid model made of hard-spheres with patchy square-well attractions. In its simplest 
version, the single patch case, the model only depends on the surface coverage x of the patch 
and the attraction range. Between the two extreme cases x = 0, the hard-sphere model, and 
X = 1, the hard-sphere square- well model, where the particles pair-potential is isotropic, the 
particles interaction is directional. The x = 1/2 model is known as the Janus the 
particle, like the roman God, has two faces of different functionalities. 

Another important process, which may lead to the manifestation of macroscopic phe- 
nomena, in certain fluids, is the clustering or association. In 1956, for example, Leon N. 
Cooper- found that the stable state of the degenerate electron fluid in a metal is one in 
which particles of opposite spin and opposite momentum form pairs. It was then under- 
stood that whereas the electrons in a metal form pairs with relative angular momentum 
zero, in ^He this would be prevented by the hard core repulsion, and that therefore Cooper 
pairing had to occur in a state of finite angular momentum. In 1961 A. Lenard^ proved an- 
alytically that a two-component plasma living in one dimension undergoes a transition from 
the conducting to the insulating state by the formation of neutral dimers made of a positive 
and a negative charge. A two-component plasma living in two dimensions is only stable at 
sufficiently high temperatures^. But if one adds a hard core to the charges it remains stable 
even at low temperatures where it undergoes the same transition^. The hard core gives rise 
to anyonic statistics for the quantum fluid living in two dimensions'^. In three dimensions 
the two-component plasma with a hard core, the so called restricted-primitive model, also 
undergoes the clustering transition at low temperature and low densities^. An example of 
a one-component Janus fluid undergoing association is the dipolar hard-sphere fluid. Here 
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a particle can be viewed as the superposition of two uniformly charged spheres: a positive 
one and a negative one^. 

In their study of the Kern and Frenkel single patch % = 1/2 Janus case, F. Sciortino 
et ai.— found that the gas branch of the coexistence curve bends at high densities at low 
temperatures. Below the critical point, the fluid tends to remain in the gas phase for a 
larger interval of densities respect to the x = 1 case. This behaviour is due to the tendency 
of particles to associate due to the directional attractive component in the pair-potential 
and form clusters. At low temperatures, these clusters interact weakly amongst themselves 
because the particles of which they are composed tend to expose the hard-sphere hemisphere 
on the outside of the collapsed cluster. 

By studying the clustering properties of the gas phase of the Janus fluid, F. Sciortino et 
ai. discovered that below the critical temperature there is a range of temperatures where 
there is formation of two kinds of preferred clusters: the micelles and the vesicles. In the 
former the particles tend to arrange themselves into a spherical shell and in the latter they 
tend to arrange themselves as two concentric spherical shells. 

It is important to confront existing cluster theories with these new findings based on com- 
puter experiments. In this work the Bjerrum cluster theory for electrolytes, later extended 
by A. Tanii^ to include trimers, has been employed (preliminary results appeared in Ref. 



15l ) for the description of the exact equilibrium cluster concentrations found in the computer 
experiment of F. Sciortino et ai.. The theory is extended to clusters of up to 12 particles in 
an attempt to reproduce the micellisation phenomenon observed in the simulations around 
a reduced temperature of 0.27. A different determination of the intra-cluster configurational 
partition function has been devised in place of the one used by J. K. Lee^. 

The Kern and Frenkel fluid has been used to describe soft matter— biological and non- 
biological materials like globular proteins in solution^'^^'^ and colloidal suspensions^'22^ or 
molecular liquids^i. Recently there has been a tremendous development in the techniques 
for the synthesis of patchy colloidal particles^^ in the laboratory. These are particles with 
dimensions of 10 — lO'^A in diameter, which obey to Boltzmann statistics^. From the realm 
of patchy colloidal particles stems the family of Janus particles for their simplicity^. It is 
possible to create Janus particles in the laboratory in large quantities'^ and to study their 
clustering properties'^. 

The micelles and the vesicles are complex structures observed in the chemistry of sur- 



factant molecules analogous to those which may be found in the physical biology of the 
cell21. 

The paper is organized as follows: in Section |Tl] we describe the fluid model, in Section 
mil we present the clustering properties of the fluid found in the Monte Carlo simulations of 
F. Sciortino et al, the cluster theory is presented and developed in Sections HVl and IVl in 
Section rvT] we compare the numerical results from our approximation to the exact results of 
F. Sciortino et al, and Section [VIII is for final remarks. 



II. THE KERN AND FRENKEL MODEL 

As in the work of F. Sciortino et al— we use the Kern and Frenkel^ single patch hard- 
sphere model of the Janus fluid. Two spherical particles attract via a square-well potential 
only if the line joining the centers of the two spheres intercepts the patch on the surfaces of 
both particles. The pair-potential is separated as follows 

$(1,2) = 0(ri2)^(ni,n2,fi2) , (1) 

where 

+00 r < cr 

-e a <r < Xa (2) 
Xa < r 

V 

and 

^ ^ , I 1 if fii ■ ri2 > cos and -112 ■ f 12 > cos 6^0 , , 

^'(ni,n2,ri2) = < (3) 
I otherwise 

where ^0 is the angular semi-amplitude of the patch. Here iii{uji) are versors pointing from 
the center of sphere i to the center of the attractive patch, with Ui their solid angles and 
ri2(fi) is the versor pointing from the center of sphere 1 to the center of sphere 2, with Q its 
solid angle. We denote with a the hard core diameter and X = 1 + A/a with A the width 
of the attractive well. 

A particle configuration is determined by its position and its orientation. 

We will use a as the unit of length and e as the unit of energy. 



One can determine the fraction of the particle surface covered by the attractive patch as 
follows 



where (. . .)a; = / • • • duj /{An). 

As in the work of F. Sciortino et ai.— we limit ourselves to the Janus case x = 1/2. 



III. CLUSTERING PROPERTIES 

The Janus fluid just described will undergo clustering as there is a directional attractive 
component in the interaction between its particles. Moreover at low temperatures the col- 
lapsed clusters are expected to interact weakly with each other. This is responsible for the 
bending at high density of the low temperature gas branch of the gas-liquid binodal curve 
recently determined in Ref. Below the critical temperature, in the vapour phase, the 
appearance of weakly interacting clusters destabilizes the liquid phase in favour of the gas 
phase. F. Sciortino et ai. during their canonical ensemble (at fixed number of particles A^, 
volume V, and temperature T, with p = N/V the density) Monte Carlo simulations of the 
fluid also studied its clustering properties. In particular they used the following topological 
definition of a cluster: an ensemble of n particles form a cluster when, starting from one 
particle, is possible to reach all other particles through a path. The path being allowed to 
move from one particle to another if there is attraction between the two particles. During 
the simulation of the fluid they counted the number Nn of clusters of n particles, which 
depends on the particles configurations, and took a statistical average of this number. 

We show in Fig. [T]the results they obtained for A = (j/2 at a reduced density pa^ = 0.01 
and various reduced temperatures ksT/e. From the figure we can see how at a reduced 
temperature of 0.27, in the vapour phase, there is the formation of two kinds of preferred 
clusters: one made up of around 10 particles and one made up of around 40 particles. 




(4) 
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FIG. 1. Exact cluster concentrations of the Janus fluid with A = (t/2 at a reduced density 
po"^ = 0.01 and various reduced temperatures ksT/e, from the Monte Carlo simulation of F. 
Sciortino et aJ.— . 

In their collapsed shape, expected at low temperatures, the particles in the clusters tend to 
expose their inactive hemisphere on the outside of the cluster, resulting in a weak interaction 
between pairs of clusters. 

In the clusters of around 10 particles the particles tend to arrange themselves into a 
spherical shell, forming a micellar structure. In the clusters of around 40 particles the 
particles are arranged into two concentric spherical shells, forming a vesicular structure. 

The aim of the present work is to see if we can approximate the exact equilibrium cluster 
concentrations found in the simulation using a cluster theory. We will restrict ourselves to 
clusters made of up to 12 particles to see if the theory is able to reproduce the micellisation 
phenomenon. The theory is described next. 

IV. A CLUSTER THEORY FOR JANUS PARTICLES 

Following Ref. [l^ we describe the fluid of particles undergoing clustering as a mixture 
of species of clusters. Clusters of species n = 1,...,N, which we call n-clusters, are 
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made up of n particles. We denote with Nn the number of clusters of species n and with 
Pn = Nn/V their density. We assume that the chemical potentials of all the cluster species 
are zero (there is no cost in energy in the formation or destruction of a cluster). Then the 
grand-canonical partition function of the fluid can be written as 

/ TV 
{N„}n=l 

where one separates the coordinates and momenta relative to the center of mass of a cluster 
from the ones of the center of mass so that q^^^"^^ will be the intra-cluster partition function 
of the cluster of species n and Qintcr the inter-cluster partition function where we consider 
the clusters as non identical. The prime indicates that the sum is restricted by the condition 
that the number of particles of the fluid is A^, 

N 

J2^^n = N. (6) 

n=l 

We approximate Qtot assuming that the sum can be replaced by its largest dominant 
contribution. Using the Stirling approximation A^! ^ {N/e)^ one then obtains 

iV 

In Qtot ~ Yl 1^ - 1^ - ^«)] + 1^ ^-t- • 

n=l 

The maximum of In Qtot as a function of {Nn} on the constraint of Eq. ([6]) is given by the 
point {Nn} where the gradients of In Qtot and of the constraint have the same direction. 
Introducing a Lagrange multiplier A the equihbrium cluster distribution {Nn} is then found 
from the conditions 

d 
dNn 

The resulting Helmholtz free energy, /3-Ftot = —In Qtot? can then be written in terms of 
the intra-cluster free energy, fif™^^^ = — Ing^"*'''^, and the inter-cluster partition function as 
follows 

N N N 

^ = E [p„ In Pn - P J + E Pn(^fn''^ + P Ju \/ - - lu Qi„ter • (9) 

n=l n=l n=l 

where /3 = l/ksT with ks Boltzmann constant and p„ = Nn/V. 

We expect the equilibrium cluster concentrations, Nn/N, to approximate the ones mea- 
sured in the simulation, {Nn)/N. 



InQt 



+ lnA" = 0, n = 1,2,3,... (8) 

{Ar„=lV„} 



V. RELATIONSHIP BETWEEN THE CONFIGURATIONAL PARTITION 
FUNCTIONS 



We will assume that Eq. also holds at the level of the configurational partition 
functions Z, as follows 

/ N 

= E n ]^ i'n''^f" Winter {{Nn} , V, T) . (10) 
{N„}n=l 

In the calculation we only work at the level of the configurational partition functions. 

Since we expect the clusters to be weakly interacting amongst themselves we will ap- 
proximate the inter-clusters configurational partition function with: i. the ideal gas approx- 
imation for pointwise clusters and ii. the Carnahan-Starling approximation^^ for clusters of 
diameter ctq- A third possibility, that we have not investigated, would be to use the Boublfk, 
Mansoori, Carnahan, and Starling approximation^i^^ for clusters of different diameters cr„. 

We will only work with a limited number ric of different cluster species. Since we are 
investigating whether the cluster theory is able to reproduce the micellisation phenomenon 
we will only consider the first Uc clusters: n = 1, 2, 3, . . . , n^. And choosing Uc = 12. 

We will describe next the two approximations used for the inter-cluster configurational 
partition function. 

A. Ideal gas approximation 

The simplest possibility is to approximate the mixture of clusters as an ideal one so that 

Zinter = ; (H) 

where Nt = Yin total number of clusters. 

The equations for the equilibrium numbers of clusters are 

N,, = X^Vz';^"^\ n= 1,2,3,..., (12) 
N = Y,nNn, (13) 

n 

from which we can determine all the concentrations Nn/N and the Lagrange multiplier 
by solving the resulting algebraic equation of order Uc. The case rzc = 2 is described in 
Appendix |X1 
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B. Carnahan-Starling approximation 



A better approximation is found if we use as the inter-cluster configurational partition 
function the Carnahan-Starling expressions^ for hard-spheres of diameter (Tq, 

In Z,„,,, = N,\riV- N,^^^^ , (14) 

where rjt = {n/Q)ptaQ is the clusters packing fraction and pt = Nt/V their density. 
In this case one needs to solve a system of ric + 1 coupled transcendental equations, 

iV„ = \^Vz^'-^^G{r)^) , z = 1, 2, 3, . . . , (15) 

N = Y,nNn. (16) 

n 

with % = {'n/Q)-ptal Pt = Nt/V, Nt = En^n, and 



G{x) = exp 



(l-.)3 ■ (^^) 

In order to search for the correct root of this system of equations it is important to 
choose the one that is continuously obtained from the physical solution of the ideal gas 
approximation as ctq — 0. Giving a volume to the clusters we introduce correlations between 
them which will prove to be essential for a qualitative reproduction of the micellisation 
phenomenon though the cluster theory. The Carnahan-Starling approximation amounts to 
choosing for the sequence of virial coefficients of the hard-spheres, a general term which 
is a particular second order polynomial and to determine the polynomial coefficients that 
approximate the third virial coefficient by its closest integer—. It could be interesting to 
repeat the calculation using for the inter-cluster partition function the hard-spheres one 
choosing all but the first virial coefficient equal to zero, to see if that is sufficient to reproduce 
the micellisation phenomena. 

Note that in order to study the vesicles we would have to solve a system of around 40 
coupled equations. 

We will describe next how do we determine the intra-cluster configurational partition 
function z^^^^^. 

C. The intra-cluster configurational partition function 

To estimate the intra-cluster configurational partition function we performed Monte Carlo 
simulations of an isolated topological cluster. 
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We determined the reduced excess internal energy per particle of the n-cluster uf^ = 
{J27<j "^(^5 = by definition) as a function of the temperature, and then used 

thermodynamic integration to determine the intra-cluster configurational partition function. 

We found that the results for u'^iT*) can be fitted by a Gaussian as follows 

<(n = a„e-^"^*' + c„, (18) 



with T* = ksT/e the reduced temperature. 

Given the excess free energy of the n-cluster iT'^^j.mtra^ then determine /^^'''^ 

l]pex,mtrs.^^ aS foUoWS 



^0 



erf (y't/^ 



(19) 



with = 1/T* and vq = ttctq/G the volume of the n-cluster. Then the intra-cluster config- 
urational partition function is given by z^^^'^'^ = exp(— n/^^'™*"^^) with z^^^'^'^ = Vq. 

We studied only the first 10 clusters with n = 3, . . . , 12. The dimer being trivial. To this 
end we started with an initial configuration of two pentagons with particles at their vertexes 
juxtaposed one above the other. The two pentagons are parallel to the {x,y) plane, have 
the z axis passing through their centers, and are placed one aX z = +cr/2 and the other at 
z = —a/2. The particles patches all point towards the origin. We formed the clusters with 
a lower number of particles by simply deleting particles and the clusters with 11 and 12 
particles by adding a particle on the z axis just above the upper pentagon and just below 
the lower one. 

We performed the simulations of the isolated cluster at a fictitious reduced density of 
pa'^ = 0.05 which ensured a simulation box big enough that the cluster did not percolate 
through the periodic boundary conditions. We also compared our results for the excess 
internal energy calculation for the isolated cluster with the results of F. Sciortino et al. for 
the low density Janus fluid, from which one extracts cluster information by taking all the 
clusters found with the same number of particles and averaging their properties, as shown 
in Fig. H 
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kgT/e 

FIG. 2. Reduced excess internal energy per particle as a function of temperature for the 6-, 9-, 
and 12-cluster. The results from the isolated (I) cluster calculation are compared with the results 
of Sciortino (S) for the Janus fluid with A = cj/2 at a reduced density pa^ = 0.01. Also shown is 
the Gaussian fit of Eq. (|18p . 

At high temperatures the limiting value for the excess internal energy per particle of the 
isolated n-cluster is —e{n — l)/n corresponding to the stretched cluster. At low temperature 
[T* < 0.15) the cluster tends to freeze into certain energy minima. So in order to reach the 
absolute minimum we used the following smoothing procedure. We smoothed the Kern and 
Frenkel potential by choosing 

^(fii,n2,ri2) = {tanh[/(ni • ri2 - cos^o)] + l}{tanh[/(-n2 ■ ri2 - cos^o)] + l}/4 . (20) 

We then gradually changed the parameter /, during the simulation, starting from 1/2 and 
increasing up to values where there is no actual difference between the smoothed potential 
and the original stepwise one. The reduced excess internal energy per particle and gyration 
radii for such minimum energy configurations are shown in Table [B 
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0.76 
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-2.50 


-15 


0.75 


7 


-2.71 


-19 


0.91 


8 


-2.88 


-23 


0.93 


9 


-3.10 


-28 


0.96 


10 


-3.20 


-32 


1.00 


11 


-3.36 


-37 


1.04 


12 


-3.42 


-41 


1.08 



TABLE I. The low temperature reduced excess internal energy per particle {U)/{en) (U is the 
potential energy of the cluster) of the clusters with up to 12 particles when A = a/2. Also shown 
is the gyration radius Rg = Y17=i ~ ^cm\'^ /n with Tcm = Y17=i being the position of the 

i-th particle in the cluster. 

In the Metropolis algorithm^i used to sample the probability distribution function pro- 
portional to e"^*^, where U is the potential energy of the cluster, the random walk moves 
through the configuration space of the particles forming the cluster through two kinds of 
moves: a displacement of the particle position and a rotation of the particle (through the 
Marsaglia algorithn>2^) . We followed two different strategies in the simulations: i. we aver- 
aged only over the particles configurations that form a cluster and ii. we explicitly modified 
the acceptance probability by rejecting moves that break the cluster. So in the second 
strategy all the moves are counted in the averages. The two strategies turned out to give 
the same results, as they should. The second strategy is preferable to simulate the bigger 
clusters at high temperature and for small well widths because there is no loss of statistics. 

In Appendix [B] we present the results for the reduced excess internal energy of the isolated 
clusters as a function of temperature and their fit of Eq. (fT8|) . 

D. Thermodynamic quantities 

Once the equilibrium cluster distribution {Nn} has been determined (within the ideal 
gas or the Carnahan-Starling approximation for the inter-cluster partition function) the 
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/3F^^ = -ln(9^) , (21) 



configurational partition function Ztot is known. Then the excess free energy is 

Z 

yN 

the reduced internal energy per particle of the fluid is 

3 1 3 ^Nndlnzf" 3 iV, 

U — ' — _ I X „ 



2(3* N d(3* 2/3* ^ d(3* 2/3* ^ iV 

n n 

and its compressibility factor, in the Carnahan-Starling approximation for the inter-cluster 
configurational partition function, is 

^ ^ 1 gin Ztot _ 1 gin Winter _ l+m + Vt~Vt (23) 

Here we have used the approximation N ^ Nt which turns out to be reasonable at the 
chosen value of the cluster diameter, as shown in Fig. HI 

In Fig. [6]we show the results for the compressibility factor and the reduced excess internal 
energy per particle. The reduced excess internal energy is compared with the Monte Carlo 
data of F. Sciortino et al. (Fig. 1 in Ref. Il3l ). 



VI. RESULTS 

We present here the numerical results from the cluster theory and compare them with the 
results of F. Sciortino et al. from the simulation of the Janus (x = 1/2) fluid with A = a/ 2. 

We studied three different attraction ranges: A = cr/2, A = cr/4, and A = O.lScr. To the 
best of our knowledge there are no Monte Carlo results available for the two smaller ranges. 

We only present the results obtained from the Carnahan-Starling approximation for the 
inter-cluster partition function as the ideal gas approximation turned out to be too crude 
an approximation even for a qualitative description of the exact clustering properties. 



A. A = cr/2 

For A = cr/2 we found the following results. 



1. Equilibrium cluster concentrations 

In Fig. [3] we compare the Monte Carlo data of F. Sciortino et al. (the results reported in 
Fig. [1]) and our results from the cluster theory. From the figure one can see that the ideal gas 
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approximation for the inter-cluster partition function is not appropriate even at high tem- 
peratures in the single fluid phase above the critical point. In order to find agreement with 
the Monte Carlo data at high temperatures it is sufficient to give a volume to the clusters, 
treating them as hard-spheres of a diameter (Tq. In the Carnahan-Starling approximation 
we gradually increased ctq from zero and found that for ctq = 2.64cr the results of the cluster 
theory were in good agreement with the Monte Carlo data at ksT /e = 0.5. Using the same 
cluster diameter at all other temperatures, we saw that the theory is able to qualitatively 
reproduce the micellisation phenomenon observed in the simulation of F. Sciortino et al. 

1 

0.1 
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0.001 
^fl 0.0001 
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2 4 6 8 10 12 

n 

FIG. 3. Comparison between the Monte Carlo (MC) data (points) and the Carnahan-StarHng (CS) 
approximation with cto = 2.64cr (lines) for the cluster concentrations {Nn)/N, n = 1, 2, 3, . . . , 12, 
as a function of the cluster size n at pa^ = 0.01 and various temperatures. Also shown is the ideal 
gas (IG) approximation at the same density and the highest temperature ksT/e = 0.5. 
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The results also suggest that with a temperature-dependent cluster diameter, or more 
generally with a cluster diameter dependent on the thermodynamic state of the fluid, we 
could achieve better agreement between our approximation and the exact results. Our topo- 
logical deflnition of a cluster has no direct geometrical interpretation. Other deflnitions with 
a geometrical nature are possible. For example Lee et al. in their studies of nucleation deflne 
an assembly of particles to be a cluster if they all lie within a sphere of radius o"o centered on 
one of the particles. In our simulations of the isolated clusters these have a globular shape 
at low temperature and a necklace shape at high temperature. The optimal cluster diameter 
ctq = 2.64(7 (found to give good agreement between the exact and approximate clusters con- 
centrations at high temperature) suggests necklace clusters made up of around 3 particles 
or globular clusters made up of around 27r(cro/c)^/"\/3 ~ 25 particles placed on a spherical 
shell. Since ctq is the only free parameter of the theory, it is important to estimate how 
thermodynamic quantities like the compressibility factor 13 P/ p, the reduced internal energy 
per particle u, and the logarithm of the total conflgurational partition function per number 
of particles. In Ztot/N, or per number of clusters. In Ztot/Nt, are sensible to variations in (Tq. 
From Fig. |4]we can see that for the thermodynamic state pa"^ = 0.01 and ksT/e = 0.5, the 
thermodynamic quantities are roughly independent of (Tq ioi ctq ^ 3cr. 
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FIG. 4. The compressibility factor, the internal energy per particle, and the logarithm of the total 
partition function per total number of particles and per total number of clusters as a function of 
the clusters diameter ctq at the thermodynamic state pa^ = 0.01 and ksT/e = 0.5 for A = O.So". 



In Fig. |5] we show the behaviour of the equilibrium cluster concentrations, from the 
Carnahan-Starling approximation with ao = 2.64cr, as a function of density at ksT jt = 0.27. 
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FIG. 5. The equilibrium cluster concentrations Nn/N, n = 1, 2, 3, . . . , 12, as a function of den- 
sity for ksT/e = 0.27 (top panel) and ksT/e = 0.5 (bottom panel) as obtained from the CS 
approximation with cjo = 2.64(7. Here A = a /2. 



From the figure we can see that at very low densities there are essentially no clusters. But 
as the density increases, clusters of an increasing number of particles appear in the fluid. 
In particular, at /cfiT/e = 0.27 there is an interval of densities where clusters of 11 particles 
are preferred. 



2. Thermodynamic quantities 

Following Section IV Dl we now use the cluster theory within the Carnahan-Starling ap- 
proximation with (To = 2.64(7 to extract thermodynamic information for the Janus fluid. In 
Fig. [6] we show the results obtained for the excess reduced internal energy per particle and 
the compressibility factor. 
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FIG. 6. The top panel shows the reduced excess internal energy per particle for three different 
values of temperature as a function of density. The results from the Carnahan-Starling (CS) 
approximation are compared with the Monte Carlo (MC) results of F. Sciortino et ai.— . The 
bottom panel shows the compressibility factor for the same values of temperature as a function of 
density from the CS approximation (no MC data is available). 



From the figure we see that there is a quahtative agreement between the results of the 
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cluster theory and the Monte Carlo results. No Monte Carlo results are available for the 
compressibility factor. 



B. A = a/4 

Decreasing the width of the attractive well to A = cr/4 yielded the results shown in Fig. 
[71 We see that now, at the reduced temperature 0.27, the preferred clusters are the ones 
made up of 10 particles. 
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FIG. 7. Same as Fig. [5] for A = cj/4. 

C. A = 0.15cr 

Decreasing the width of the attractive well even further to A = O.lScr, we obtained the 
results of Fig. |8l Now, at the reduced temperature 0.27, there is a range of densities around 
pa^ = 0.1 where the preferred clusters are made up of 7 or 8 particles. 
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FIG. 8. Same as Fig. [5] for A = O.lSa. 

VII. CONCLUSIONS 

We constructed a cluster theory for a fluid undergoing clustering and showed that it 
is able to reproduce the micellisation phenomena recently observed in the simulation of 
the vapour phase of Kern and Frenkel Janus particles^. A topological definition of the 
cluster is used. We determined the intra-cluster configurational partition function through 
thermodynamic integration of the excess internal energy of the cluster, estimated through 
Monte Carlo simulations of an isolated cluster. In the simulation we restricted the random 
walk through the configurations of the particles that compose the cluster by rejecting the 
moves that break the cluster. Due to the geometrical characteristics of the pair-potential it 
is expected that the clusters, when in their collapsed shape, will be very weakly interacting 
amongst themselves as the Janus particles will expose the hard-sphere hemisphere on the 
outside of the cluster. We thus used for the estimation of the inter-cluster configurational 
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partition function first the simple ideal gas approximation for pointwise clusters and then 
the Carnahan-Starling approximation for clusters seen as hard-spheres of diameter (Tq. The 
equilibrium cluster concentrations obtained with the ideal gas approximation turned out 
to disagree, even at high temperatures, with the ones obtained from the simulation of the 
fluid^ and were not able to reproduce the micellisation phenomenon in the vapour phase. 
We then gradually increased ctq from zero until we found good agreement between the 
equilibrium cluster concentrations obtained with the Carnahan-Starling approximation and 
the concentrations from the simulation of the fluid^^ at high temperature (above the critical 
point). Using the same value of (Tq for lower temperatures (below the critical point) we were 
able to qualitatively reproduce the micellisation phenomenon observed in the simulation of 
the fluidi^ around a reduced temperature of 0.27 and a reduced density of 0.01. This result 
is important for two reasons. Firstly it shows that the clustering fundamentally arising 
from the canonical ensemble description of the fluid of particles can be approximated by a 
grand canonical ensemble description of a particular clustered fluid. Secondly the second 
description, which assumes from the start a clustered structure of the fluid, is much less 
computationally costly than the first. Unlike most previous works on cluster theories where 
the aim is usually to avoid the Monte Carlo simulation^i^, our approach is a hybrid one 
where we still use the Monte Carlo experiment to determine the intra-cluster properties. Of 
course our goal can only be a qualitative description of the fluid as we specifically prescribe 
a particular description of the clusters and this is the source of our approximation. 

Studying the behaviour of the equilibrium cluster concentrations as a function of den- 
sity and temperature, we saw that the micellisation phenomenon only takes place within 
a particular range of temperatures (below the critical point) and densities (in the vapour 
phase) . 

Once the equilibrium concentrations have been found it is possible to determine how the 
cluster theory approximates the thermodynamic quantities of the fluid. We find qualitative 
agreement between the Monte Carlo data of F. Sciortino et ah— and our approximation for 
the excess internal energy of the vapour phase. For the compressibility factor no Monte 
Carlo data is available so our results remain a theoretical prediction. 

We studied three different values of the attractive square- well width: A = a/2, A = a/A, 
and A = 0.15a. Monte Carlo results^ are available only for the largest width. Our study 
shows that as the range of the attraction diminishes the micelles tend to be made up of a 
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smaller number of particles. 

A related interesting problem to that just discussed is the one of trying to give a definition 
of a liquid drop expected to form in the coexistence region as a result of the condensation 
instability. 



Appendix A: Connection with Wertheim association theory 

At small Xi allowing only clusters of one (monomers) and two (dimers) particles, we get 

iVi = AV^i"*"^" , (Al) 
iVs = X^Vzf'" , (A2) 
N = Ni + 2N2 , (A3) 

which is a quadratic equation in A. The solution for the fraction of patches that are not 
bonded (fraction of monomers) is 

^ = A=^. (A4) 

P 1 + Vl + 8pA 

with A = /[z'^^'^'^f and p = N/V the density of the fluid, in accord, at low T, with the 
recent analysis of F. Sciortino et al— (compare their X of Eq. (10) with our Pi/p and their 
A with our A), based on Wertheim association theory^^. Our theory, contrary to the one of 
Wertheim, allows to consider the case of multiple bonding of the patch. 



At high temperature our A differs from the A of Ref. 
begin to dissociate. 



34 but in this limit the clusters 



Appendix B: Tables for the excess internal energy per particle of the clusters 

We present here the results for the reduced excess internal energy per particle as a function 
of temperature of the isolated n-cluster with n = 2, 3, . . . , 12 as obtained from our Monte 
Carlo simulations. 

In Tab. Ullwe show the results at A = 0.5a obtained with the strategy i. described in Sec. 
IV C[ The smoothing procedure described in Sec. IV CI was used at the lowest temperature. 
The excess internal energy per particle of the n = 2 cluster is always — e/2 given our 
topological definition of a cluster. 
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In Tab. IIIII we show the results at A = 0.5a obtained with strategy ii. described in Sec. 
IV CI The smoothing procedure described in Sec. IV CI was not used at the lowest temperature. 
Comparing Tabs. Hi] and HIT] we can see that the two strategies lead to the same results. 

In Tab. IIVI we show the results at A = 0.25cr obtained with strategy i. described in Sec. 
IV CI The smoothing procedure described in Sec. IV CI was not used at the lowest temperature. 

In Tab. |V]we show the results at A = O.lScr obtained with strategy ii. described in Sec. 
IV CI The smoothing procedure described in Sec. IV CI was not used at the lowest temperature. 

In Tab. I VI I we give the fit to the Gaussian of Eq. ( ITSj) of the reduced excess internal 
energy per particle as a function of the temperature. 
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TABLE II. The tables refer from left to right to clusters made up of n = 3, 4, . . . , 12 particles. U 



is the potential energy of a cluster of n particles. Below ksT/e = 0.1 the reduced excess internal 
energy per particle remains roughly constant in all cases: the smoothing procedure described in 
Sec. IV CI was used. The data was obtained with a Monte Carlo simulation over 5 million steps 
where one step consists of n particles moves. The strategy i. described in Sec. IV CI was used in 
the simulations. 
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TABLE III. The tables refer, from left to right, to clusters made up of n = 3,4, 10, 11 particles. U 
is the potential energy of a cluster of n particles. The smoothing procedure described in Sec. IV CI 
was not used at the lowest temperature. The strategy ii. described in Sec. IV CI was used in the 
simulations. 





<(7>/(En) 


error 




k 


bT/ 


e (U) / (en) error 


ksT/e {U)/{tn) 


error 






{U)/{en 


error 




k 


BT/e {U)/{in) 


error 


oo 


-0.666 









oo 


-0.75 





oo 


-0.8 







oo 


-0.833 









oo 


-0.857 





0.7 


-0.705 


0.002 






0.7 








0.7 


-0.87 


0.02 




0.7 














0.7 






0.5 


-0.732 


0.002 






0.5 


-0.866 0.007 


0.5 


-1.00 


0.03 




0.5 


-0.95 


0.01 






0.5 


-0.95 


0.01 


0.3 


-0.832 


0.002 






0.3 


-1.138 0.005 


0.3 


-1.427 


0.008 




0.3 


-1.63 


0.01 






0.3 


-1.79 


0.01 


0.1 


-0.99872 


0.00008 






0.1 


-1.4987 0.0002 


0.1 


-1.7984 


0.0002 




0.1 


-2.16 


56 


0.0002 






0.1 


-2.22 


0.02 


kBT/t 


(U)/(en) 


error 






e (U) / (^n) error 




feflT/E {U)/{en 


error 






E ((7)/{«l) 


error 




fe 


BT/e ((7)/{«i 


error 






oo 


-0.875 







oo 




0.888 




oo 


-0.9 







oo 




0.909 









oo 


0.916 









0.7 








0.7 








0.7 








0.7 












0.7 












0.5 








0.5 








0.5 








0.5 












0.5 












0.3 


-1.91 


0.03 




0.3 




-1.95 0.02 




0.3 


-2.07 


0.04 




0.3 




-2.10 


0.04 






0.3 


2.01 


0.03 






0.1 


-2.3706 


0.0009 




0.1 




2.4416 0.0005 




0.1 


-2.5969 


0.0006 




0.1 




2.721 


0.002 






0.1 


2.730 


0.008 







TABLE IV. Same as Table [IT] but with A = 0.25fj. The smoothing procedure described in Sec. 



IV CI was not used at the lowest temperature. 
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TABLE V. Same as Table HIl but with A = 0.15cj. The smoothing procedure described in Sec. IV CI 
was not used at the lowest temperature. The strategy ii. described in Sec. IV CI was used in the 
simulation. 
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TABLE VL Fit to the Gaussian of Eq. (jlSp of the reduced excess internal energy per particle of 
the first eleven n-clusters as a function of temperature. 
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